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ABSTRACT 

The timescale for energy release is an important parameter for constraining 
the coronal heating mechanism. Observations of "warm" coronal loops (~ 1 MK) 
have indicated that the heating is impulsive and that coronal plasma is far from 
equilibrium. In contrast, observations at higher temperatures (~ 3MK) have 
generally been consistent with steady heating models. Previous observations, 
however, have not been able to exclude the possibility that the high temperature 
loops are actually composed of many small scale threads that are in various stages 
of heating and cooling and only appear to be in equilibrium. With new obser- 
vations from the EUV Imaging Spectrometer (EIS) and X-ray Telescope (XRT) 
on Hinode we have the ability to investigate the properties of high temperature 
coronal plasma in extraordinary detail. We examine the emission in the core of 
an active region and find three independent lines of evidence for steady heating. 
We find that the emission observed in XRT is generally steady for hours, with a 
fluctuation level of approximately 15% in an individual pixel. Short-lived impul- 
sive heating events are observed, but they appear to be unrelated to the steady 
emission that dominates the active region. Furthermore, we find no evidence for 
warm emission that is spatially correlated with the hot emission, as would be ex- 
pected if the high temperature loops are the result of impulsive heating. Finally, 
we also find that intensities in the "moss," the footpoints of high temperature 
loops, are consistent with steady heating models provided that we account for the 
local expansion of the loop from the base of the transition region to the corona. 
In combination, these results provide strong evidence that the heating in the core 
of an active region is effectively steady, that is, the time between heating events 
is short relative to the relevant radiative and conductive cooling times. 
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Introduction 



The coronal heating problem is one of the most fundamental open questions in solar 
physics. The formation of the high temperature solar atmosphere is clearly related to the 
magnetic fields generated in the solar interior, but how this magnetic energy is converted 
to the thermal energy of the corona remains unknown. One important constraint on the 
coronal heating mechanism is the time scale for energy release. If energy is deposited into 
magnetic flux tubes on timescales that are very short compared to a characteristic cooling 
time, then the corona will be filled with loops that are close to equilibrium and appear to 
be steady. Active region observations ha ve generally suggested tha t high temperature loop s 
are consistent with steady heating (e.g., IPorter fc Klimchukl Il995l : iKano fc Tsunet all 19951 ). 
These studies have found that the observed evolution of the emission is much slower than the 
radiative and conductive cooling timescales. T here has also been some success in modeling 



entire active regions with steady heating models dSchriiver et al 
20061 : Iwinebarger et al.ll2008l : iLundquist et al.ll2008h . Finally, 



2004; 



Warren & Wine 



Antiochos et al 



barger 



(|2003|) have 



argued that observations of the "moss," the footpoints of high temperature loops, are also 
consistent with steady heating. They found that the average moss intensities are typically 
constant over many hours and loops cooling through 1 MK were not observed in the moss 
region they observed. 

There is considerable evidence that coronal loops observed at lowe r temperatures ( 
1 MK) are evolving and not in equilibrium (e 



20031 : ICirtain et all 120071 : lUgarte-Urra et al 



200 



Aschwanden et al 

IT 



2001 



Winebarger et al 



These "warm" loops appear to have 
apex densities that are much higher than can be acco unted for by steady heating models 
(e.g. , IWinebarger et al.l 120031 : I Aschwanden et al.l 120081 ) . The properties of the warm loops 



are more consistent with impulsive heating models (e.g., IWarren et al J 12002 ; ISpadaro et al. 
20031 : IWarren et al.ll2003h . 



If it is true that hot loops are close to equilibrium while the warm loops are generally 
heated impulsively then the coronal heating mechanism becomes even more difficult to un- 
derstand. The cooling time for short, hot loops is relatively rapid (typically a few hundred 
seconds), so heating events on these loops would need to occur very frequently. The cooling 
time for the warm loops is much longer (typically a few thousand seconds), so impulsive 
heating events on these loops would be infrequent. 

It is tempting to conjecture that the emission at high temperatures is also consistent 
with impulsive heating models and that the apparent steadiness of this emission is the result 
of th e supe r position of many evolving stand s along the line of sight (e.g. ICargill fc Klimchuk 
19971 , I2004J : IPatsourakos fc Klimchukl 120061 ) . Previous observations have not been able to 
exclude this possibility. For example, it could be that the cooling loops are not easy to 
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detect in an active region core because they are faint relative to both the bright moss and 
the extended corona in which the active region is embedded. 

The EUV Imaging Spectrometer (EIS) and the X-ray Telescope (XRT) on the Hinode 
mission provide a new opportunity to observe active region emission in unprecedented detail. 
EIS observes emission covering a very broad range of coronal temperatures: between Fe VII 
and Fe XXIV, only 4 ionization stages of Fe are not present in the EIS data (Fe XVIII-xxi). 
EIS also observes Ca XIV, XV, xvi, and xvii, pro viding excellent coverage of the critical 
temperature range around 3 MK ( IWarren et al XRT is a broadband imaging telescope 



that observes high temperature plasma very efficiently and at high spatial resolution. The 
high cadence XRT data complement EIS, which often observes at a much lower cadence, and 
allows us to track the evolution of coronal plasma over a large field of view. 

In this paper we use EIS and XRT observations to examine the properties of coronal 
plasma in the core of an active region. The active region that we have selected (NOAA active 
region 10960) is unusual in that most of the overlying warm loops are located to the north 
and south of the region, providing a largely unobstructed view of the active region core over 
a wide range of temperatures (see Figured]). High cadence EIS observations of this region, 
which was observed by Hinode during the period 4-13 June 2007, have shown that the Fe XII 
195.119 A intensities, Doppler shifts, and non-th ermal widths in the mos s are constant over 



long periods of time, suggesting steady heating ( IBrooks fe Warrenl 12009 ) . In this study we 
find three additional lines of evidence that also indicate that the heating of high temperature 
loops is steady. We examine the evolution of the emission in individual pixels in XRT and 
find that the vast majority of the emission is constant to within approximately 15% during 
this period. The inspection of emission at different temperatures in the core of the active 
region shows that there is no evidence for loops cooling from high temperatures. We find 
no relationship between the warm emission (Fe X-Fe XIV) and the steady hot emission 
(Ca XlV-Ca XVII). Finally, we find that we can bring all of the observed moss intensities 
into agreement with steady heating models, if we allow for loop constriction at the base of the 
loop. Previous mod eling work had foun d significant discrepancies at the lowest temperatures 



observed with EIS ( IWarren et al.ll2008l ). 



None of these observations is conclusive; they do not necessarily exclude alternative 
models that involve non-equilibrium processes. However, these observations do provide very 
strong constraints on the mechanism responsible for producing the high temperature emission 
observed in solar active regions. 
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2. Observations 

The EIS instrument on Hinode is a high spatial and spectral resolution imaging spec- 
trograph. EIS observes two wavelength ranges, 171 212 A and 245 291 A, with a spectral 
resolution of about 22 mA and a spatial resolution of about 1" per pixel. There are 1" and 
2" slits as well as 40" and 266" slots available. The slit-slot mechanism and the CCD both 
image an area on the Sun 1024" in height, but a maximum of 512 pixels on the CCD can be 
read out at one time. Solar images can be made using one of the slots or by stepping one 
of the slits over a region of the Sun. Te lemetry constraints g enera lly limit the spat i al and 



spectral coverage of an observation. See ICulhane et al.l ( 120071 ) and iKorendyke et al.l (120061 ) 
for more details on the EIS instrument. 

For these observations the 1" slit was stepped over the central part of the active region 
and 25 s exposures were taken at each position. An area of 128" x 128" was imaged in about 
57 minutes (see Figured]). The observing sequence for this observation returned the complete 
wavelength range from each detector so all of the lines observable with EIS are potentially 
available in these data. 

The raw data were processed using eis_prep to remove dark current, warm pixels, and 
other instrumental effects using standard software. During the processing the observed count 
rates are converted to physical units. Intensities from the processed data are computed by 
fitting the observed line profiles with Gaussians. One line of particular interest that requires 
special attention is Ca XVII 192.858 A, which is strongly blended with Fe XI 192.813 A 
and an O V multiplet near 192.90 A. An algorithm for deconvolving this blend using Fe XI 
188.216 A and multi-component fitting has been developed and is discussed in detail by 



Ko et al.l (120091 ). We use this approach here to determine the Ca XVII intensities. The 
uncertainties for this line are significantly greater than for the other lines that we consider. 
However, as illustrated by the EIS rasters shown in Figure Ej, there is a strong similarity 
between the deconvolved Ca XVII 192.858 A image and other rasters from ions formed at 
similar temperatures, such as Ca XV 200.972 A and Fe XVII 254.87 A. 

The XRT on Hinode is a high cadence, high spatial resolution (approximately 1" pixels) 
grazing incidence telescope that images the Sun in the soft X-ray and extreme ultraviolet 
wavelength ranges. Temperature discrimination is achieved through the use of focal plane 
filters. Because XRT can observe the Sun at short wavelengths, XRT images can observe high 
temperature solar plasma very efficiently. The thinner XRT filters allow longer wavelength 
EUV emission to be imaged and extend t he X RT response to lower temperatures. Further 



details on XRT are given in iGolub et al.l ( 120071 ) . 



The XRT data taken around the time of the EIS raster consisted of relatively high 
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cadence Open/Ti-Poly images with Open/Al-thick, Open/Be-thick, and G band images 
interleaved at a cadence of about 600 s. For this analysis we have processed all of the 
Open/Ti-Poly images with xrt_prep to remove the dark current and to do exposure time 
normalization. All of the images have been co-aligned to the initial frame using cross corre- 
lation. To investigate the longer term evolution of this region we have considered all of the 
XRT data within ±2 hours of the start of the EIS raster. An example XRT image is given in 
Figured! An animation of these data is available in the electronic version of the manuscript. 

In this analysis we also m ake reference to ob servations from TRACE \ the Transition 



Region and Coronal Explorer ( IHandy et al.lll999l ). TRACE is a normal incidence, multi- 



layer telescope. There are 3 channels for imaging the corona: FelX/x 171 A, FeXII 195 A, 
and Fe XV 284 A. There is also a long wavelength channel for imaging the photosphere, 
chromosphere, and transition region. All of the images are projected onto a single CCD and 
images in different wavelengths must be taken sequentially. For these observations the data 
are mostly from the 171 A channel, with a few white light, 1600 A, and 284 A images taken 
for context. The 171 A data were processed with trace_prep, despiked, and co-aligned using 
cross-correlation. An example TRACE 171 A image is given in Figure [TJ An animation of 
the data is available in the electronic version of the manuscript. 



3. Analysis 

In this section we presented a detailed analysis of the EIS, XRT, and TRACE data 
available for this period. The EIS data represent a snapshot of the active region accumu- 
lated over about an hour and yield little insight into the temporal evolution of the emission. 
The XRT and TRACE observations, in contrast, have excellent temporal coverage but some- 
what limited diagnostic capabilities. In combination they allow us to establish a detailed 
understanding of the plasma properties in the core of an active region and how they evolve 
in time. 



3.1. Temporal Evolution 

As is suggested by the four-hour XRT movie associated with Figure EH the large scale 
pattern of the soft X-ray emission from this region is remarkably steady. To illustrate this 
we have selected several points from the core of the active region and plotted light curves 
for the emission in these individual pixels. The light curves, shown in Figure El generally 
show fluctuation levels of about 10-15% around the median intensity for this time interval. 
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In Figure [3] we also show two points that illustrate the evolution of transient loops. The 
lifetimes of these brightenings are generally on the order of 1000 s or less and show that the 
intensities in the core of the active region are steady on timescales that are long compared 
to a cooling time. 

A more systematic view of the variability is obtained by computing the median intensity 
(I) and standard deviation (07) for each point in the co-aligned XRT data cube. A spatially 
resolved plot of oj/ J is shown in Figure [H Note we have considered only those pixels with 
median intensities above 50 DN s _1 . Low intensities well outside the core of the active region 
are highly variable. 

This calculation shows that there a few structures in the core of the active region that 
are highly variable, such as those illustrated in Figure El bu t that the em ission in most 
of the active region is constant to within approximately 15%. IShimizul ( 119951 ) performed a 
comprehensive study of active region transient brightenings with the Soft X-ray Telescope on 
Yohkoh. He found that the distribution of event energies followed a power-law distribution 
and that the total energy in these events was about a factor of 5 smaller than the energy 
required to heat the active corona. Our result is qualitatively consistent with this systematic 
study. 

The XRT light curves suggest that the emission in the core of an active region is 
largely steady, with relatively few significant brightenings and a very high median inten- 
sity level (> 1000 DN s" 1 ). This result is consistent with previous results from SXT (e.g., 
Porter fc Klimchukl 1 19951 : iKano fc Tsunetal 1 19951 ) , which showed that the observed emission 
decayed on timescales that were long relative to a cooling time. To put these observations 
in perspective, the cooling time for a 50 Mm loop cooling from equilibrium at 5MK is of the 
order of 800 s, while the XRT emission is relatively steady for several hours. Note that a 
50 Mm loop would connect the middle of the moss regions in this active region. 

The core of the active region is clearly made up of the emission from many different loops 
and it is unclear if the low fluctuation levels are the result of steady heating; on ind i vidual 
loops or from the superpositio n of many evolving loops ( jCargill fc Klimchukl Il997l . 12004] : 
Patsourakos fc Klimchukl 120061 ) . Models based on many evolving loops, however, predict 
that there should be warm emission (~ 1 MK) t hat is spatially correlated wit h the cooling of 
hot loops in the core of the active region (e.g., IWarren fc Winebargerl 120071 ) . Furthermore, 
the magnitude of the warm emission should be related to the observed intensity of the hot 
emission. Rising levels of hot emission should be accompanied by a similar increase in the 
warm emission. These predictions can be tested easily with the observations of this region 
with EIS. 
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3.2. Warm and Hot Emission 

The interpretation of many active region observations is complicated by the presence 
of overlying warm loops. For this region such loops are largely, but not completely, absent, 
giving a relatively unobstructed view of the core of the active region at warm temperatures. 
This is illustrated by the TRACE 171 A movies associated with Figure [TJ In these movies 
we see large scale 1 MK loops to the north and south of the active region core but not over 
the core itself. We also do not see many loops cooling through the TRACE 171 A bandpass 



that a re connected to the moss, similar to the observations considered by lAntiochos et al. 



(|2003|). 



The limited amount of overlying emission allows us to look for any faint warm emission 
that is related to the formation of hot coronal loops. Another useful aspect of these obser- 
vations is that there is a strong gradient in the intensities of the hot lines (Ca XIV-xvn, 
Fe XVII, and XRT) in the core of the active region in between the moss. To illustrate the 
behavior of the emission at different temperatures in the core of the active region we have 
selected a line segment parallel to the moss and extracted the intensities from each of the 
EIS rasters. These intensities are displayed in Figure [5] where we compare the intensities in 
each line with those from a hot line, Ca XV 200.972 A. To facilitate the comparisons we have 
normalized the intensities to their values at the end of the line segment where the intensities 
of the hot lines are generally the smallest. For the hot lines the intensity along the segment 
increases by a factor of 3 to 4 in this region. At the cooler temperatures (Fe X-xiv), however, 
the intensities are unchanged and show no spatial correlation with the increase of the hot 
emission. The behavior of the emission from Fe XV and Fe XVI appears to fall in between 
the hot and warm lines, with Fe XVI following Ca XV fairly closely and Fe XV behaving 
more like a warm line. 

Since the EIS raster lacks temporal information it is possible to argue that the more 
intense emission in the hot lines is evolving and will show up at a later time as enhanced 
emission at warmer temperatures. We have, however, already used the XRT observations to 
demonstrate that hot emission in the core is constant over many hours (Figure [3]). Inspection 
of later XRT data for this active region shows that the strong gr adients in the intensity 



and the dark inter-moss region actually persist for several days (see Brooks fe Warren 



2009 



Figures 1 and 3). The hot emission shown in the EIS rasters is clearly not dominated by 
transient events. 

We note that the intensities of the warm lines are non-zero in the core of the active 
region. Inspection of the TRACE image shown in Figure [1] suggests that while the inter- 
moss region is the dimmest part of the active region core, it is not the dimmest part of the 
active region at these temperatures. As we will show in the next section, almost all of the 
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warm lines have inter-moss intensities that are less than 20% of the moss intensities. The 
Fe XIII 202.044 A line, however, is unusually bright in the inter-moss region. The lower level 
of the transition responsible for this line is very sensitive to collisonal de-excitation into other 
levels and the emissivity for this line rises sharply with decreasing density. This sensitivity 
to low density background plasma may account for the enhanced intensity of the Fe XIII 
202.044 A line in the inter-moss region. 

The origin of the warm emission in the core of the active region is unclear. It is 
most likely related to unresolved, high lying loops that form around the active region (e.g. 
Mason et al.lll999l ). Scattered light from the bright moss may also influence the observed in- 
tensities in this region. See, for example, iDeForest et al.l ( 120091 ) for a discussion of scattered 
light in the TRACE instrument. Scattered light levels have not yet been measured in EIS, 
but because EIS and TRACE have several similar design elements, such as a front entrance 
filter supported by a mesh and multi-layer coatings, observations of dim regions with EIS 
are also likely to suffer from some level of stray light. 



4. Modeling the Moss 



The moss represents the footpoints of high temperature loops and offers important 
boundary conditions for physical models. Observations of the moss are particularly useful 
for constraining physical models since the ob served intensity is p r oportional to the loo p 
pressure and independent of the loop length ( Martens et al.l l2000l : IVourlidas et al.l l200ll ) . 
That is, 

h oc P x / (1) 

where I\ is the observed intensity, P is the base pressure, and / is the filling factor. This 
allows for the calculation of physical models of the moss without knowing the loop length. 
This is important since the magnetic topol ogy of an active re gion is difficult to infer, even 
with the use of vector magnetograms (e.g.. LDeRosa et al.ll2009l ). 



Some previous work comparing EIS moss intensitie s with steady heating models has 



been presented in an earlier paper ( 1 Warren et al.l 120081 ) . This work introduced the use 
of density sensitive lines to determine both the base pressure and the filling factor. For 
lines formed close to Fe XII there was good agreement between steady heating models and 
the observations. At the lowest temperatures, however, there was a dramatic difference, 
with the modeled intensities being about 400% higher than what was observed. A similar 
temperature dependen c e in t he contrast between moss and network intensities was noted by 
Fletcher fe de Pontieu Their work suggested that it was difficult to identify moss 

regions in relatively low temperature emission. 
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Here we investigate the possibility that variations in the loop cross section are responsible 
for the relatively low moss intensities at low temperatures. For some time it has been 
recognized that magnetic flux tubes must expand rapidly in the region between the high 
beta photosphere and the low beta corona. That is, coronal loops must look like "funnels" 
in the transition region. This loop exp ansion plays an important role in determinin g the 



r 

energy balance on coronal loops (e.g., iGabriell 119761 : iDowdy et al.l 119871 : iRabinl 119911) and 



needs to be accounted for in modeling the moss intensities. The network model of [Gabriel 
( 119761 ) shows a significant constriction at temperatures below about 1 MK, qualitatively 
similar to what has been suggested by the EIS observations of the moss. 

The first step in performing detailed comparisons between theory and observation with 
these data is to accu mulate the obse r ved m oss intensities from this active region. As shown 
in Figure E2 we follow I Warren et al.l ( 120081 ) and use a density sensitive line ratio to identify 
the dense moss regions. For this work we focus on the Fe XIII 203.826/202.044 A ratio. The 
threshold has been set at logn e = 9.7 so that density contours encompass the bright region 
in the core of the active region. For each line we have extracted the intensities in this region 
and computed the distribution of intensities. The resulting distribution is approximately 
Gaussian and for each line we have determined the median intensity as well as the standard 
deviation in the distribution of intensities. A summary of all of the observed moss intensities 
is given in Table [TJ 

To account for any contamination from overlying loops or scattered light, we have sub- 
tracted an approximate background derived from the low intensity pixels in the inter-moss 
region. For consistency, the same background subtraction procedure is applied to all of the 
emission lines. The background intensities are given in Table [TJ 

It is unclear what uncertainties should be associated with these median intensities. The 
moss intensities are relatively high and the statistical errors in the fits to the Gaussian line 
profiles are generally small. Furthermore, for simplicity we will work with the median moss 
intensities and the uncertainties in these values are also small. Finally, the errors in our 
analysis are undoubtedly dominated by the systematic errors in the atomic data and by the 
assumptions embedded in the hydro modeling, such the assumption of ionization equilibrium. 
Absent any compelling alternative we assume that the uncertainty in each measured intensity 
is 15%, which is much higher than the statistical uncertainty. 

Note that in a previous paper we used the Fe XII 186.880/195.119 lin e pair to identify the 
moss and determine the pressure and filling factor (jWarren et al.ll2008l ). Recent work, how- 
ever, has shown tha t there are discrepa ncies between the various Fe XII and Fe XIII density 
sensitive line ratios (jYoung et al.ll2009l ). leading one to wonder which densities are the most 
accurate. As we will discuss in more detail in the next section, the Si X 258.375/261.058 A 
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ratio is generally consistent with the results from Fe XIII and we use the Fe XIII ratio here. 

If we assume a constant loop cross section and steady, uniform heating, the only re- 
maining parameters needed to solve the hydrodynamic loop equations are the volumetric 
heating rate and the loop length. As indicated by Equation [Q any family of solutions with 
the same base pressure will yield the same observed moss intensity. To illustrate this we cal- 
culate a grid of solutions that covers a range of loop lengths and heating rates. To solve the 
hydrostatic loop equations we use a numerical code written by Aad van Ballegooijen (e.g., 



Schrijver fc van Ballegooijenl 120051 ). We consider total loop lengths in the range L = 10- 
100 Mm and heating rates that yield maximum temperatures in the range \ogT max = 2.5- 
7.5 MK. In the numerical model the lower boundary condition is set at 0.02 MK and the 
loops are assumed to be oriented perpendicular to the solar surface. 

For each numerical solution we calculate the total intensity in the line using the usual 
expression 

Ix = ^LJ 6A ( Te ' ne ) n e ds ^ ( 2 ) 

where 6\{ T P , n P ) is the emissivity computed fr om the CHIANTI atomic physics database ver- 
sion 5.2.1 jDere et al.l[l997l : lLandi et alJbood ). Figure H shows the calculated 203.826/202.044 A 



ratio as a function of the pressure at 1.5 MK, the temperature of formation for Fe XIII. The 
ratio is clearly a function of the pressure alone and allows us to easily convert the observed 
line ratio into a value for the base pressure. An example calculation for the Fe XI 188.216 A 
intensity is also shown in Figure [3 With the pressure determined we can now read off the 
calculated intensities for each line. For Fe XI 188.216 A, for example, the calculated intensity 
is 8625 erg cm -2 s _1 sr _1 , which is much larger than what is observed and indicates that the 
emission is not resolved at the spatial resolution of EIS. 

Calculated intensities for all of the emission lines are given in Table [TJ We have estimated 
a filling factor of 20% using the median ratio of the observed to modeled intensities for the 
Fe X, XI, XII, and xill lines. This corresponds to the Fe XI 188.216 A line. For the emission 
lines from these ions there is generally good agreement between the observed and modeled 
intensities. At lower temperatures, however, the observed intensities are much lower than 
what is calculated from the simulation. Ratios of the observed to modeled intensities as a 
function of temperature of formation are displayed in Figure El This plot shows that the 
constant area model breaks down at temperatures be low about 1 MK. T his discrepancy is 



consistent with the comparisons presented previously ( IWarren et al.l 120081 ). 



An alternative approach to modeling the observed emission is to assume a loop length 
and determine the heating rate and filling factor that best fits the observed intensities. This 
method is based on the assumption that the observed moss intensities are independent of the 
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loop length, an assumption supported by theory (IMartens et al.ll2000l ) and by the numerical 
solutions presented here. This method has the advantage of being extensible. Our goal is 
to consider loops with complex geometries and the intensity-pressure approach outlined in 
the previous paragraphs becomes very cumbersome when applied to models with several free 
parameters. We have implemented such an approach using the Levenberg-Marquardt least- 
squares minimization routine MPFIT for the Interactive Data Language (IDL). The resulting 
modeled intensities are presented in Table [1] and are very similar to those derived from the 
intensity-pressure relationships. We have assumed a loop length of 50 Mm and obtain best-fit 
parameters of eo = 5.98 x 10 _3 erg cm -3 s _1 and / = 0.22. Only the emission lines formed 
at Fe X and above are considered in the minimization. In this approach all of the lines 
of interest are considered simultaneously and the pressure indicated by the Fe XIII lines is 
not necessarily reproduced by the model. To emphasize these lines in the final solution we 
reduce the uncertainties assumed for these lines to 5%. 

To explore the role of geometry in determining the observed intensities we parametrize 
the cross-sectional area as a function of height using the following expression 



1 + 



r- 1 



tanh 



s- s 



+ 1 



(3) 



for s < L/2. The area is extended to the other footpoint by symmetry. The area expansion 
is defined so that the minimum is A(s — 0) ~ 1/T and the maximum is A(s = L/2) ~ 1. 
The gradient in area expansion peaks at Sq and the magnitude of the gradient is inversely 
proportional to a s . Smaller values of a s indicate a sharper transition. This form has been 
assumed because it allows for the loop expansion to be localized. Previous work has consid- 
ered functions of the form A ~ T a , where the expansion continues throughout the corona 



(e.g., lRabinlll99ll : IChae et al.lll998l ). T here is consider able empirical evidence that coronal 
loops have constant cross sections (e.g.. iKlimchukl 120001 ) and our area function is consistent 
with this. Note that in our prescription the loop geometry is fixed and does not respond to 
heating within the loop. 

The expectation is that the constriction will reduce the calculated intensities at the low 
temperatures. The discrepancy between the observed intensities appears to occur around the 
temperature of formation of Fe X, which is approximately 1 MK. In the constant cross-section 
case this temperatures occurs at about 0.5 Mm up the loop length. 

To solve for the best-fit parameters we assume a loop length of 50 Mm and use the 
Levenberg-Marquardt approach. The free parameters are the volumetric heating rate (e ), 
the filling factor (/), and the loop expansion parameters T, s , and a s . An example solution 
is shown in Figure [9] and the corresponding calculated intensities are given in Table [TJ 
The solutions to the constant cross section and funnel models are generally similar. The 
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constant cross section model has somewhat higher temperatures and lower densities. The 
most significant impact on the calculated intensities comes from the area term in the emission 
measure 

EM = A(s)n 2 e ds. (4) 

The funnel model clearly reproduces the intensities at the lower temperatures much better 
than the constant cross-section models do. In this fit all of the emission lines are used. To 
emphasize the Fe XIII lines over the other lines the error assumed for the Fe XIII intensities 
is again reduced to 5%. 



Rabin 


1991; 


Chae et al. 


1998) 



have considered the area as a function of temperature and for comparisons with these works 
we show A(T(s)) derived from the model parameters and the solution to the loop equations 
in Figure flOl O ur result is qualitatively similar to the most promising funnel model derived 
by Rabin (fl 99 ll . see model "B4HL" in their Figure 12). In this model the area is relatively 
constant until logT ~ 5.6 and then expands rapidly, similar to wha t we show in Fig ure [9l 



Our result, however, is not consistent with the model derived by IChae et al.l ( 119981 ) from 
SUMER Doppler shift measurements (see their Figure 19). Their model shows significant 
expansion at temperatures above logT ~ 5.0. In both cases, as the authors were well aware, 
comparisons were made with spatially averaged quiet Sun observations a nd the applicability 



of the se results to our modeling of the active region moss is unclear. IPatsourakos et al. 



(1 19991 ) used spatially resolved measurements of the quiet network to examine loop expansion 
from the transition region to the corona. The amount of loop expansion that they measured 
is generally similar to what we have determined here, although it appears to begin at a 
somewhat lower temperature than in our model. 

One factor that w e have not consid e red is the impact of chromospheric absorption on 
the moss intensities. De Pontieu et al. (1999) have shown that the evolving, reticulated 
pattern of the moss is often related to obscuration by chromospheric jets on neighboring 
field lines. Absorption leads to reduced observed intensities and impacts the filling factors 
and heating rates determined from the observations. This chromospheric absorption could 
also be stronger at the lowest heights and could explain the lower intensities of the lower 
temperature lines. This could also potentially complicate the inference of the transition 
region geometry from the observed intensities. Recent analysis of coordinated EIS and 
SUMER observatio ns suggests that emiss ion below the Lyman continuum is reduced by 
about a factor of 2 JDe Pontieu et al.ll2009h . It should be noted that this result is predicated 
on comparisons between 2 Fe XII line ratios (186.880/195.119 A and 1241.990/195.119 A). As 
we have discussed, the atomic data for Fe XII is problematic and there is some uncertainty 
in this correction. The results from Fe XII also do not yield insights into the possibility of 
additional absorption at lower heights. Clearly more analysis is needed on this important 
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issue. 



5. Electron Densities 

In this paper we have used the Fe XIII 203.826/202.044 A line ratio to determine the 
electron density while previously the Fe XII 186.880/195.119 A ratio was used. The Fe XII 
ratio is sensitive over a very wide range of densities (logn e = 7 - 12), making it a very 
useful plasma diagnost ic. It has been fou nd, however, that the Fe XIII and Fe XII densities 



generally don't agree ( jYoung et al.ll2009h . leading to the potential for systematic errors in 
the modeling. 

In many previous studies the atomi c data for Si has proven to be very robust (e.g, 



Feldman et al.l Il999l : iDoschek et al.l 119971 ) and so we have compared the densities derived 
from the Si X 258.375/261.058 A ratio with those from the other line pairs. The Si X ratio 
is sensitive over a much smaller range of densities than the Fe XII and Fe XIII pairs so we 
have looked for observations of all three line pairs in small active regions, bright points, and 
the quiet Sun, i.e., regions where the densities are generally lower than in the moss region 
that we have studied here. 

The result of one such comparison is shown in Figure [TT1 This plot shows the densities 
computed in a small active region from all three line pairs. These calculations clearly show 
that at high densities the Fe XII results are systematically higher than those from Fe XIII 
and Si X and that the Fe XIII and Si X densities are in good agreement. We have found 
similar results from the analysis of 4 other quiet Sun and active region observations. 

The problems with Fe XII become more pronounced at the highest densities, which is 
particularly problematic for the analysis of the moss. In these observations the peak density 
is about a factor of 3 higher in Fe XII. Since the observed intensity scales with the square of 
the density this discrepancy is amplified. We have derived filling factors of 10-20% for this 
moss region. Emphasizing the Fe XII ratio would lead to filling factors almost an order of 
magnitude smaller. 



6. Summary 

We have presented a comprehensive analysis of observations in the core of an active 
region using data from the EIS and XRT instruments on Hinode and TRACE. The apparent 
steadiness of the XRT emission, the lack of spatial correlation between the hot and warm 
emission, and the consistency of the funnel models with the observed emission all point to 
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frequent heating events that keep the hot loops close to equilibrium. Furthermore, these 
results are consistent with high cadence EIS measurements of moss intensities, Doppler 
shifts, and nonthermal width s that show little evidence of dynamical events over many 



hours (IBrooks fc Warrenl 120091 ) . In combination, these results provide strong evidence that 
the heating in the core of an active region is effectively steady, that is, the time between 
heating events is short relative to the relevant radiative and conductive cooling times. 

Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ 
as domestic partner and NASA and STFC (UK) as international partners. It is operated 
by these agencies in co-operation with ESA and NSC (Norway). TRACE is supported by 
a NASA contract to LMATC. The authors would like to thank Jim Klimchuk for helpful 
discussions on the coronal heating problem, George Doschek for an explanation of the Fe XIII 
202.044 A emissivity, and the referee for a number of very helpful comments on the original 
manuscript. 
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Fig. 1. — Context images for observations of NOAA active region 10960 taken on 5 June 
2007. Top left panel: An MDl/SoHO magnetogram with the TRACE and XRT fields of 
view indicated. Top right and bottom left panels: TRACE 171 A and XRT Ti-Poly images 
with the EIS field of view indicated. Bottom right panel: A composite TRACE and XRT 
image showing the difference in morphology between the hot and warm emission. Animations 
of the TRACE and XRT data are available with the electronic version of the manuscript 
(moviel_trace.mpg, movie l_trace_moss.mpg and movie2_xrt.mpg). 
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Fig. 2. — EIS rasters from the core of NOAA active region 10960 taken on 5 Jun 2007 between 
06:48:20 and 07:44:59 UT. With the exception of Ca XVII 192.858 A, which was deconvolved 
from a blend with other lines using a procedure described in the text, the intensities were 
derived from simple single or double Gaussian fits to the line profile. 
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Fig. 3. — The evolution of the high temperature emission measured at various points in the 
active region core. The left panels show a single image from the XRT data set with the 
points indicated. The right panels show the light curves for these points. For each point 
the median intensity (7) is indicated with a solid line and / ± aj is indicated by the dashed 
lines. The red bar at the top of panel highlights the time of the EIS raster. The vertical 
line indicates the time of the XRT image displayed on the left. The top panel illustrates 
relatively constant emission that dominates the core of the active region. The bottom panels 
illustrate the evolution of more variable emission. 
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Fig. 4. — Plots of I and a//7 as a function of position for the XRT Open/Ti-Poly data. For 
the oij I plot only the points with a median intensity above 50 DN s _1 are shown. For low 
intensity pixels the fluctuation level is generally very high (> 30%). 
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Fig. 5. — Intensities in the core of an active region. Top panels: EIS rasters in Fe X, Fe XII, 
Fe XVI, Ca XIV, and Ca XVII. Bottom panels: The intensities along a segment in various 
emission lines. For each plot the intensities are compared with Ca XV 200.972 A. At the 
lowest temperatures (Fe X-Fe XV) there is little or no relationship between the hot and 
warm emission. At the highest temperatures (Fe XVI-Ca XVII) the intensities are strongly 
correlated. The intensities are normalized to the values at the end of the segment (~ 40"). 
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Fig. 6. — The moss identified using a simple threshold with the density derived from the 
Fe XIII 203.826/202.044 A ratio. The contour is for log n e = 9.7. The box indicates the region 
used to estimate the background contribution to the line intensities. For each emission line 
considered here we have computed the median intensity in the moss and background regions. 




Fig. 7 — Top panel: The Fe XIII 203.826/202.044 A ratio as a function of the pressure at 
1.5 MK calculated from a steady heating model. Solutions over a range of loop lengths and 
volumetric heating rates have been used. The red dot represents the pressure corresponding 
to the median intensities given in Table [U Bottom Panel: The intensity in the Fe XI 
188.216 A line as a function of pressure. 
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Fig. 8. — The ratio of observed to calculated intensity for the expanding area (or funnel) 
model and constant area models (const2) plotted as a function of temperature. This plot 
illustrates the failure of the constant cross-section models at low temperatures. 



- 25 - 




Fig. 9. — The transition region loop geometry derived from EIS observations of the Moss. 
Top panel: The normalized area as a function of distance along the loop. The best-fit 
parameters are e = 8.03 x 10~ 3 , / = 0.11, T = 45.6, s = 1.1 Mm, and cr s = 1.6 Mm. 
Bottom panels: The temperature and density along the loop. The peak temperature of 
formation for each ion is indicated on the plots. For comparison, the corresponding best-fit 
solution for the constant cross section case is also plotted. The total loop length is 50 Mm 
and only part of the solution is shown. 
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Fig. 10.— The area aS 3j function of temperature for the expanding loop model shown in 
Figure M The peak temperature of formation for each ion is indicated. 
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Fig. 11. — EIS observations of a small active region beginning on 2 December 2007 at 
11:10:44 UT. Top panels: Rasters in 3 pairs of density sensitive lines. Intensities have been 
extracted from the region indicated by the white line. Bottom panels: Densities computed 
from the 3 pairs of lines. The final panel shows the densities as a function of position. The 
Fe XIII and Si X densities are in generally good agreement while the Fe XII densities are 
systematically higher. 
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Table 1. Observed and Modeled Moss Intensities 8, 



Ion 


-^median 


cri 




Ibac 


k 


-^observed 


Iconstl 


R 


Iconst2 


R 


1 funnel 


R 


Mg V 276.579 


36.4 


12, 


.4 


3. 


,1 


33.3 


178, 


.3 


0. 


.19 


231 


.3 


0. 


.14 


38. 


,0 





.88 


Mg VI 270.394 


97.8 


25 


.9 


8. 


,5 


89.3 


420. 


.1 


0, 


.21 


541, 


.3 


0. 


.16 


95. 


A 





.94 


Mg VII 280.737 


96.7 


27, 


.5 


5. 


,1 


91.6 


247, 


.1 





.37 


321 


.4 


0. 


.28 


67. 


A 


1 


.36 


Si VII 275.352 


136.6 


42, 


.0 


9. 


,6 


126.9 


529 


.6 


0. 


.24 


668. 


.2 


0. 


.19 


137, 


.7 





.92 


Fe IX 188.497 


248.0 


64, 


.6 


13. 


,8 


234.2 


710, 


.1 





.33 


892, 


.6 


0, 


.26 


229. 


.9 


1 


.02 


Fe IX 197.862 


106.1 


26 


.1 


8. 


,0 


98.0 


310. 


.3 


0. 


.32 


384. 


.8 


0, 


.25 


101, 


.8 





.96 


Fe X 184.536 


936.9 


225 


.3 


62. 


.3 


874.6 


1107. 


.0 





.79 


1352 


.1 


0. 


.65 


553, 


.3 


1 


.58 


FeXI 182.167 


1114.6 


266, 


.3 


59. 


.0 


1055.6 


1114, 


.1 


0, 


.95 


1326 


.6 


0. 


.80 


935. 


.0 


1 


.13 


Fe XI 188.216 


1827.7 


398, 


.5 


121, 


.4 


1706.4 


1706, 


A 


1 


.00 


1984 


.4 


0. 


.86 


1347. 


.6 


1 


.27 


Si X 258.375 


1390.2 


301, 


.4 


65. 


,4 


1324.8 


1171. 


.6 


1 


.13 


1354 


.4 


0. 


.98 


1355. 


.8 





.98 


Si X 261.058 


455.0 


97. 


.3 


24. 


.7 


430.3 


355 


.0 


1 


.21 


406, 


.9 


1, 


.06 


402. 


.0 


1 


.07 


Fe XII 186.880 


3026.4 


676, 


.5 


130. 


.8 


2895.6 


2608. 


.4 


1 


.11 


3072 


.6 


0. 


.94 


3263, 


.0 





.89 


Fe XII 192.394 


1309.5 


266, 


.7 


81. 


,9 


1227.5 


1443. 


.1 





.85 


1622 


.4 


0. 


.76 


1571, 


.6 





.78 


FeXII 195.119 


4102.9 


827, 


.4 


271. 


.6 


3831.3 


4883, 


.6 





.78 


5516 


.1 


0. 


.69 


5395. 


.9 





.71 


Fe XIII 202.044 


2196.5 


349, 


.7 


214. 


,8 


1981.7 


1457, 


.3 


1 


.36 


1626 


.5 


1, 


.22 


1759. 


.2 


1 


.13 


Fe XIII 203.826 


6609.6 


1284 


.1 


304, 


,9 


6304.7 


4633 


.7 


1 


.36 


5379 


.5 


1, 


.17 


6446. 


.1 





.98 



a In this table wavelengths are in A and the intensities are in erg cm -2 s _1 sr _1 . Note that 07 is the Gaussian 
width for the distribution of intensities in the moss and not the statistical uncertainty, hack is an estimate of 
the background emission in the moss. I CO nsti is the computed intensity for the constant cross-section model 
computed from the intensity-pressure relationships, I CO nst2 is the constant cross-section model computed from 
the least-squares fitting, and I funnel is the best-fit model for the funnel geometry. R is the ratio of observed to 
calculated intensities. 



